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Abstract: In this paper, a new method is proposed to automatically detect screw threads in 3D 
density fields obtained from computed tomography measurement devices. The described method 
can be used to automate many operations during screw thread inspection process and drastically 
reduce operator's influence on the measurement process resulting in lower measurement times 
and increased repeatability 
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1. INTRODUCTION 



Many types of coordinate measurement systems are widely 
used to inspect parameters of metric screw threads (NPL, 
2010). Screw thread measurement using coordinate mea- 
suring machines (CMM) is a tedious process which requires 
high operator's skill. CMM measurements usually aquire a 
limited set of points suitable for a feature-based inspection, 
however not enough for an in-depth analysis and assesment 
of a screw pair fit. This limitation is crucial while dealing 
with screw threads for mission critical applications (petro- 
chemical industry, heavy engineering industry et al.). Now- 
days, as a result of well-established international practice, 
the complex inspection of screw threads is performed using 
thread gauges. Gauges are subject to wear and tear and 
require regular replacement and inspection which yields 
expenses. 

In this paper, a new method is proposed, allowing to 
perform feature-based inspection and complex in-depth 
analysis of screw threads using the date aquired from 
computed tomography (CT) measuring devices. Experi- 
mental part of this work was done using X-ray CT scanner 
Metrotom 1500 (see Fig.l) by Carl Zeiss IMT GmbH. 

2. RELATED WORK 

Threaded and plain gauges are widely used for screw 
thread inspection during recent decades. They allow com- 
plex inspection of threaded fits (NPL, 2010). Gauges are 
expensive inspection instrument since they are subject to 
wear and the whole inspection procedure can be long- 
lasting for large diameter threads (M150 and above). 
Thread gauges that are proved to be out-of-range by direct 
measurements can be fit with good master gauges. This 
situation is typical but gauges and direct measurements 
are not mutually exclusive methods. Considering the NPL 
experience (NPL, 2010) one can use both master gauges 

* Presented images were rendered using Linderdaum Engine, 
http://www.linderdaum.com. 



and direct measurements to do an elementwise inspection 
of the thread gauges. First, one should measure pitch, 
pitch diameter and inner/outer diameters of the thread. 
Thereupon gauges should be applied to assist the direct 
measurements. Elementwise inspection is more accurate 
and can yield numerical results (instead of simple go/no- 
go answer). The point is direct elementwise measurements 
do not evaluate the whole surface of the thread therefore 
form deviations can be overlooked. Gauges can prevent 
this kind of errors even after elementwise measurements 
succeeded. 

Gauges are unable to separate pitch and diameter devia- 
tions. In addition, pitch errors can be hidden by increased 
pitch diameter. It is traditionally considered practical 
to use gauges for small threads. As described in NPL 
notes (NPL, 2010) the primary reason for this discrepancy 
between gauges and direct measurements is thread form 
deviations. 

Nowadays, a lot of elementwise thread assessment methods 
have evolved. Their majority is based on coordinate mea- 
surement devices and numeric evaluation of the results. 
A lot of researchers work torwards the improvement of 
thread inspection. In (Ikonomov et al., 1995) the methos 
of position and form measurement is proposed based vir- 
tual gauges. Carmignato and Chiffre (Carmignato and 
De Chiffre, 2003) proposed a screw inspection method with 
a special needle-like probe fitted on a coordinate measuring 
machine. In (He et al., 2006) authors perform optical 
inspection of damaged screw threads using CCD camera. 
Many techniques usually involve different shape analysis 
algorithms from the field of computer vision and pattern 
recognition. Robertson and Fisher (Robertson and Fisher, 
2001) experimented with 3D scanners and their applica- 
tion to large thread measurement. They examined that it 
is possible to extract parameters of screw threads from 3D 
scanned data. However they deal only with inner and outer 
radii of the thread. In the proceedings (Katz et al., 2009) a 
method is presented to perform an elementwise inspection 
of internal threads using laser sensor and CCD camera. 
Kosarevsky and Latypov (Kosarevsky and Latypov, 2010) 




Fig. 1. X-ray CT scanner Carl Zeiss Metrotom 1500 

used Hough transformation to extract features from planar 
sections obtained via profile measuring machines. The 
mathematics behind these image recognition techniques 
can be found in (Princen et al., 1992). 



Metrotom 800 X-ray CT scanners can achieve values of 
MPE E = 4.5 + O.OLL /am (Benninger, 2009). 

3. INITIAL APPROACH 



In (Perng et al., 2010) a working system is proposed 
capable of automatic internal thread inspection. It is based 
on industrial endoscope and computer vision algorithms. 
However, the main purpose of this system is to detect 
surface defects and not its geometrical properties. An 
in-depth overview of classical thread inspection meth- 
ods is provided in NPL Notes on Screw Gauges (NPL, 
2010). Current element-wise thread inspection methods 
are mainly based on the work of G.Berndt (Berndt, 1940). 

Application of high resolution computer tomography (CT) 
in dimensional metrology has grown popular during the 
recent years. It has moved from qualitative assessment of 
workpieces in non-distructive tests to a precise measure- 
ment instrument (Kruth et al., 2011). 

The accuracy of spatial measurements in these tasks is 
highly dependent on the geometry and material work- 
pieces. A lot of numeric compensation algorithms are used 
to achieve the highest possible accuracy and reduce differ- 
ent artifacts of CT scanning. 

Traceability of coordinate measurements obtained from 
computed tomography devices is ensured via special cali- 
bration gadgets and procedures. Nowadays these methods 
received ISO certification. Modern software can extract 
geometrical features from measured density fields and 
evaluate their parameters according to ISO norms. To 
assess the accuracy of CT systems (besides calibration) 
sets of reference workpieces are used. These workpieces 
are calibrated on high-precision CMMs and results are 
compared to the CT data. Recent experiments on GE 
Sensing & Inspection Technologies (Germany) CT devices 
show (Luebbehuesen, 2009) the deviation of results for 
distances and diameters to be within 6 jam. Carl Zeiss 



In this paper, a method of planar sections is used to 
numericaly assess the quality of screw threads. It allows 
to reduce the problem of 3D shapes recognition to the 
finite number of 2D recognition problems. Our method 
operates on a bundle of planes that contain thread axis. 
The principal symmetry axis of the CT-scan is considered 
to be the thread axis. Performing image recognition in each 
plane a numeric evaluation of thread parameters can be 
obrained. The planar evaluation algorithm was presented 
in (Kosarevsky and Latypov, 2010). Using the obtained 
results one can assess the quality of the measured thread. 

Let f(xi,X2,X3) : R 3 i->> R be the material density 
distribution inside the object measured by a CT device and 
let V C R 3 be the domain where function / is defined. Let 
us consider a bundle of planes P a that share a common 
point r = (^1,^2,^3) ("mass center" of the measured 
object) as "distinguished" planar sections of the thread. 
Components r are evaluated as follows 

_ Sv x if ( x i, x 2, xs) dxi dx 2 dx 3 



Jx> f( x i-> x 2->xz) dx\ dx2 dxs 



0,1,2. (1) 



The thread axis is the common axis of planes P a which 
is the principal symmetry axis of the scanned object. 
Symmetry axes are approximated with the eigenvectors 
of the covariance matrix 



/M200 M110 M101 
A = I /ino M020 M011 
\/iioi M011 M002 . 



(2) 



where 



l^ijk 



= I w ijk 



f(xi,x 2 ,xs) dxi dx 2 dx 3 . 



W ijk = (sci - xxf (x 2 - x 2 ) J (x 3 - x 3 ) . 



(3) 




Fig. 2. Planar section passing through the axis of M5 thread 

The plane P a passing through the axis of the internal M5 
thread is shown in Fig. 2 with the section of the object. 



4. ANALYTIC REPRESENTATION OF A SCREW 
THREAD 

In (Nicolson and Fearing, 1991) and (Nicolson and Fearing, 
1993) Nicolson et al. proposed a model to represent a basic 
profile of the metric screw surface. Their model is based 
on piecewise functions and is effective at finding contact 
surfaces. In this work a method more suitable for thread 
recognition is used. 

Consider the family of helicoidal surfaces 



A bundle of coordinate surfaces intersecting the plane 
x = is shown in Fig. 4. 

Basic profile (ISO, 1998) of the metric screw thread (Fig. 5) 
is defined by the pitch P, the height of generating triangle 
H and the nominal diameter D. 

According to (ISO, 1998) the values are: 
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For a metric screw thread the profile angle is a = 60°, so 
that tan f = ^ an( ^ equations (7) yield: 
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paramerized by u E [— 1; 1], v € [0; V] and w G [0; +oo). 

Assumption w = D yields the surface in Fig. 3: 

The transform (x\ y\ z) — »■ (u\ v\ w) can be used to check if 
some point (x; y; z) belongs to the surface w = D. From (4) 



it follows that v — z. Also, x 2 + 1 
and x/y = tan (Rv + iru), which gives 
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It is easy to see that planar sections of the surface w - 
are piecewise linear. Projection on Oyz plane yields x 
that is sin (Rv + iru) = or 



(5) 
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and u is selected from the [— 1; 1] segment by the appro- 
priate choice of the k value. 
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Any longitudinal planar section of the nominal screw 
thread is the planar section of the surface Fig. 4, truncated 
on both sides according to the parameter c. 

5. ASSESMENT OF SCREW THREAD 
PARAMETERS 

The overall workflow of our algorithm is as follows: 

(1) Acquire scalar 3D density field. Authors used 
Carl Zeiss Metrotom 1500 X-Ray CT scanner to 
acquire data. Also, several sets of generated data were 
used to perform an in-depth data analyzis and proove 
the correctness of our implementation. 

(2) Perform feature extraction. Apply the threshold- 
ing filter. Non-zero values correspond to in-material 
points. Fig. 6 shows the thresholding result for the 
scanned external M5 thread with a threshold value of 
0.48. The actual range of interest is specified manu- 
ally. 

(3) Establish thread axis. Calculate covariance matrix 
using the equation (2) and determine its eigenvectors. 




Fig. 3. Surface w = D 
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Fig. 4. Section of the surface w = D at x = 



One of the eigen vectors defines as approximate (6) 
thread axis, which is a common line for all of the 
2D cross section planes. 

(4) Extract isosufrace. Apply the three-dimensional 
Sobel filter to determine the points of thread surface. 
The result of 3D Sobel operator is shown in Fig. 7. 

(5) Generate planar sections. Fix an arbitrary plane 
containing the center of mass of the volume and 
previously determined thread axis. Rotate this plane 
using a number of uniformly distributed angles in the 

. . . 27r range. Clusterize points into buckets corre- (7) 
sponsing to the closest plane. Every single point can 
get into a bucket only once. 



Process each planar section (bucket) as de- 
scribed in (Kosarevsky and Latypov, 2010). 

Straight lines representing flanks of the thread are 
obtained via Hough transformation. Sort obtained 
lines and find intersection points. Intersection points 
closest to the profile average line are used to fit circles 
into screw thread groves. Fit circles into the groves. 
Radii of the circles are chosen so that circles contact 
the profile close to the pitch diameter, as in the three- 
wires method. 

Evaluate thread pitch and pitch diameter. Dis- 
tances between the fitted circles are used for assess- 
ment of screw thread parameters in every planar 
section. 




Fig. 5. Basic profile of ISO metric screw thread (ISO, 1998) 




Fig. 6. M5 screw thread density field after the thresholding filter application 




Fig. 7. Results of 3D Sobel operator application to the region of density field 

The described procedure can be applied not only to segmentation is already done and one can continue from 
density fields, but also to 3D point clouds. In that case step 3, using the coordinates of separate points. 



Table 1. Numeric results 



Thread 


Pitch 


Pitch std. dev., 




(min.max.), 


mm 




mm 




External thread M5 


0.97 -M.01 


0.02 


Internal thread M5 


0.91-M.07 


0.06 


Artificial model 


1.00 -M.00 


^0.00 



6. EXPERIMENTS AND AUTOGENERATED 
MODELS 

The correctness of the algorithm is checked by processing 
a number of generated volumes corresponding to some 
screw threads with specified parameters. The generated 
volume data is an 1024 3 8-bit three-dimensional array of 
samples which is essentialy a quantized graph of the xs 
function corresponding to the screw thread. To estimate 
the precision of the algorithm the generated volume data 
is corrupted with gaussian noise. 

The sample of the generated screw thread model is pre- 
sented in Fig. 8. 

Two aluminum samples with metric screw threads were 
scanned using Carl Zeiss Metrotom 1500 CT X-ray scan- 
ner: internal and external M5 ISO screw threads. Numeric 
results of screw thread pitch evaluation using 36 planar 
sections are presented in the table 1. 

Significant value of standard deviation for internal thread 
is due to thicker material and size of the workpiece 
compared to the external thread workpiece. Standard 
deviation for our artificial model was below 0.001 mm. 
However, probability distribution and dependence on the 
input data need to be assessed more carefully. 

Our software was implemented in C++ using Min- 
pack (Min, 1999) optimization library and KHT Sand- 
box (https://sourceforge.net/projects/khtsandbox) (Fer- 
nandes and Oliveira, 2008), the reference implementation 
of the kernel-based Hough transform for detecting straight 
lines in binary images. It allows a software implementation 
to achieve real-time performance even on relatively large 
images. Fast evaluation of f(xi, #2, #3) gradients was done 
in a GLSL shader. 

Direct evaluation of moments /i^ using the formula (3) 
(necessary to establish thead axis), can take significant 
time at higher resolution. However, in this paper, the 
problem is not addressed. In case where material density 
corresponds to a metric screw thead 1 special optimiza- 
tions methods can be used (Flusser and Suk, 1993). 

Eigenvectors of a covariance matrix are calculated using 
Jacobi eigenvalue algorithm (Flannery et al., 2002). 

Presented images were rendered using Linderdaum Engine 
and Linderdaum Volume Rendering Library. 

7. FUTURE WORK 

Future work must be undertaken to assess the measure- 
ment uncertainty of the proposed algorithm and how 
noises in the input data and their filtering affect the 



1 Segmented volume contains only binary and 1 without any 
intermediate values. 



results. Measurement uncertainty can be determined via 
Monte-Carlo trials according to (GUM, 1995) and (GUM, 

2008). 

It is important to mention that finding values of param- 
eters P and D so that the virtual gauge will fit into the 
measured thread without intersection will enable direct 
evaluation of virtual pitch diameter of the thread. 

8. CONCLUSIONS 

One of the main advantages of the proposed method is 
that the result is not just a binary "go/not-go", but 
contains additional numeric data which can be evaluated 
to determine the exact thread features with deviations. It 
is important for further technological decisions concerning 
the production of threaded parts. 

One of the main flaws of the proposed methods is the 
scanning speed of CT equipment. It is required up to 
one hour of machine time to inspect a single workpiece. 
However, other types of surface scanning devices (i.e. 
laser scanning of external threads) can yield improved 
performance and provide usable 3D point clouds. 
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